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ABSTRACT 

We apply the K-correction to the black hole LMXB GX 339-4 which implies M x > 6M 
by only assuming that the companion is more massive than ~ 0.17Mq, the lower limit al- 
lowed by applying a 'stripped-giant' model. This evolutionary model successfully reproduces 
the observed properties of the system. We obtain a maximum mass for the companion of 
M2 < l.lM© and an upper limit to the mass ratio of q(= M^/Mx) < 0.125. The high 
X-ray activity displayed by the source suggests a relatively large mass transfer rate which, 
according to the model, results in M2 > O.3Af and Mx > 7M & . We have also applied this 
scenario to the black hole binary XTE J 1550-564, which has a similar orbital period but the 
donor is detected spectroscopically. The model successfully reproduces the observed stellar 
parameters. 
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1 INTRODUCTION 

Low Mass X-ray Binaries (LMXB) are interacting binaries har- 
bouring a neutron star (NS) or a black-hole (BH) which accretes 
matter from a low mass star. A subgroup, the so-called X-ray tran- 
sients (SXTs), provide a unique opportunity to set dynamical con- 
straints to the the masses of both NSs and BHs because the spec- 
trum of the companion dominates during the quiescent (X-ray off) 
states. In particular, we can empirically determine the mass func- 
tion of the system which is usually expressed as: 



f(M) 



2-kG 



M x sin i 3 

(i + g) 2 



(1) 



where P is the orbital period, K2 the semi-amplitude of the 
velocity curve of companion star, Mx the mass of the compact 
object, i is the orbital inclination of the systems and q — j^ 2 - the 
mass ratio. Therefore, it is possible to establish a secure lower limit 
to the mass of the compact object by only measuring P and Ki. 
Although this meth od has success fully been applied to ~ 20 BH 
X-ray binaries (see ICasaresll2006l for details), they only represent 
the 'peak of the ice-berg' of a large populatio n of 10 8 — 10 9 
stella r-mass BH present in our galaxy (e. g. iBrown & Bethel 
1994 ). Most quiescent SXTs are optically faint, with R > 20, 
and hence dynamical constraints are usually affected by large 
errors. Setting new dynamical constraints and refining the existing 
mass solutions are the only ways to get new insights into the 
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fundamental properties of stellar-mass BHs. 
GX 339-4 (V821 Ara) was discovered by OSO 7 in 1972 
dMarkert et al.lll9"73h . Since then, several X-ray outburst (e. g. 4 in 
the last decade) have been observed , and all the X-ray states have 
been detected in this system (See iMendez & van derKlij[l997l 
for the intermediate sta te andlMivamoto et alj 1991 for the Very 
High state). Moreover, ICorbel et alj j2000h detected a compact 
jet emission from this system, showing that it is a microquasar. 
GX 339-4 was ea rly proposed as a BH candidate based on its 
X-ray properties dSamimi et alj 1 19791) but the spectral features 
of the comp anion star could not be detected even during 'X-ra y 
off states dShahbaz. Fender & Charles] 12001: hereafter SFCO lb . 
preventing a dynamical confi rmation (i. e . Mx > 3Mq). Only 
during the 2002 outburst iHvnes et al] J2003; hereafter H03I) 
reported the first detection of the donor star thanks to the discovery 
of NIII/CIII Bowen emission lines arising from the irradiated 
companion star. The lines are very sharp and swing with a velocity 
semi-amplitude of K ern =3l7 ± 10 km s" 1 . These authors also 
reported an orbital per iod of P = 1.7557 ± .0004 days which 
was later confirmed bv lLevine & Corbel d2006h . The combination 
of both, radial velocity of the companion and orbital period yields 
f{M) = 5.8 ± O.5M (and hence M x > 5.3M ), which 
represents the first dynamical proof for a BH in GX 339-4. Note 
that narrow high-excitation emission lines originating from the 
compa nion star have als o been detected in many others LMXBs 
(e.g. ICasares et alj 120041) and demonstrate that this is a common 
feature in LMXBs with high X-ray activity (i. e. steady systems and 
transients during outbursts). However, the NIII/CIII emission lines 
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are excited on the inner hemisphe re of the donor star by the Bowen 
mech anism/photoionization dMcClintock. Canizares & Tarter] 
1 19751) and only provide s a lower limit (K em ) to the true K2- 
velocity of the donor. InlMufioz-Darias. Casares & Martmez-Paisl 
d2005; hereafter MCM05I) we tackle this problem in a general 
approach by modeling the deviation between the reprocessed 
light-center and the center of mass of a Roche lobe filling star (the 
so-called 'K-correction') including screening effects by a flared 
accretion disc. In this paper we compute the K-correction for GX 
339-4 which provides a more restrictive lower limit to the mass of 
the BH in this binary. In sections [3] and [4] we discuss the nature of 
the companion star and show that all the known observables can 
be explained by an scenario in which a 'stripped-giant' companion 
transfers mass onto a BH. We also apply this evolutionary model to 
XTE J1550-564 which shows observational properties very similar 
to GX 339-4. The stellar parameters of the companion in this 
system are in excellent agreement with a stripped-giant scenario. 



2 K-CORRECTION FOR GX 339-4 

Following MCM05 the relation between the radial velocity corre- 
sponding to the center of mass of the companion (K2) and the ob- 
served emission line velocity (if em ) can be expressed as: 

K 

1 < in 



K 2 = 



1-/(1+9) 



(2) 



where / is a dimensionless factor which represents the displace- 
ment of the emission light-center from the center of mass of the 
companion in units of the orbital separation. The value of / de- 
pends on how much area of the irradiated hemisphere of the donor 
is shadowed by the flared accretion disc. According to MCM05, / 
is constrained between 



0.5 + 0.227 log q > f > sin a M 



(3) 



where the left hand side of the equation represents the case when 
the emission line is formed in the inner Lagrange point whereas the 
right hand side corresponds to the case where the emission arises 
from the irradiated point with maximum radial velocity. This point 
is located at the limb of the irradiated region and holds an opening 
angle <xm above the plane of the binary. An analytical expression 
of dM as function of q is provided by Paczynski's equation (1971), 



sinaM 3* — 0.462( — 
a 1 
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(4) 



It is clear from eq. [2]that K2 (and hence f(M)) increases with /. 
From the lower limit to / in eq.|5] it is possible to establish a more 
restrictive lower limit to Mx than the one provided by H03 f(M) 
(i. e. Mx > 5.3Mq). In fig.[T]we show the lower limit obtained 
for Mx as function of q (thick solid line). Furthermore, GX 339- 
4 does not show eclipses of the X-ray source which results in the 
following restriction to the orbital inclination (i): 

i<90-a A/ (g). 

By using these constraints in eq.Q]we obtain a strict lower limit to 
Mx as function of q. This is shown in fig. Q] as a thick solid line. 
The dotted line shows the lower limit to f(M) reported by H03 
but including also the above constraint on i and the dependence of 
f(M) with (1 + q) 2 . Hence, the difference between the solid and 
dotted lines gives the net effect of the K-correction in the lower 
limit of the BH mass. For instance, it rises the lower limit by ~ 
0.5M Q for q = 0.05 and ~ 1M© for q = 0.1. Unfortunately 
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Figure 1. K-correction for GX 339-4. We show a strict lower limit to 
Mx (thick solid line) as function of the mass ratio (q) and the mass 
of the companion(A/2). We rind Mx > 6.0Mq (thick solid line) and 
q < 0.125. The dotted line shows the same limit by considering only the 
value reported by H03 (Mx > 5.3Mq), the absence of X-ray eclipses 
and the dependence of f(M) with (1 + q) 2 (see eq. [TJ. We also mark 
f(M) = 5.8 ± 0.5M a reported by H03 as the light, grey region. 



there are no strong constraints on q and hence a restriction to the 
companion mass are clearly needed in order to derive a lower limit 
to Mx independently of q. 



3 THE NATURE OF THE COMPANION STAR IN GX 
339-4 

In terms of LMXB evolution there are three possible cases in which 
stable mass transfer occurs. These three evolutive paths essen- 
tially depend on what time scale (angular momentum loss or stel- 
lar ev olution) brings the companion into Roche lobe contact first 
(e. g. iKing. Kolb,& Burderilll996; hereafter KKB9rj) . In this pic- 
ture, LMXBs with orbital periods around ~ 2 days descend from 
systems where the companion has evolved off the main sequence 
as a subgiant before angular mome ntum losses shrink the Roche 
lobe e nough to allow mass overflow. iFaulkner. Flannerv. & Warned 
d 19721) showed that the mean density (p) of a Roche lobe-filling 
companion star is determined solely by the binary period P: 

9= 113P,- 2 gcm- 3 

where P^r is the orbital period in hours. For GX 339-4, with 
P hr = 42.14 hours, we obtain p = 0.06 g cm~ s which according 
to ICoxl d200d) would correspond to a B main sequence star with 
~ 17M©. This scenario does not fit with the GX 339-4 case since 
a B star would dominate the optical spectrum of the binary. On 
the other hand, late-type giants have even lower mean densities 
than the companion in GX 339-4. Therefore, the comp anion star is 

robably a subgiant, a s it was previously suggested by Hvn es et all 

2004; hereafter H04 . 

A subgiant is a low mass star which has evolved off the main 
sequence and posses a helium core which represe nts an ap proxi- 
mate fraction ~ 0.15 of the total stellar mass (see King 1988 for 
details). rTaamldl983h and lWebbink. Rappaport. & Savonije dl983l) 
have shown that all the stellar properties depend on the core mass 
(M c ) instead of its total mass (M2). In this scenario mass transfer 
is driven by the expansion of the Hydrogen burning shell as con- 
sequence of the nuclear evolution. This so-called 'stripped-giant' 



Table 1. Stripped giant model for GX339-4 





Minimum Mass 
solution 


Maximum Mass 
solution 


Core Mass (M ) 


0.166 


0.187 


Total Mass (M ) 


0.166 


1.1 


Radius (Rq) 


1.56 


2.93 


Luminosity (Lq) 


1.19 


3.25 




4837 


4533 


-M 2 (M Q yr" 1 ) 


4.9 x 10~ n 


7.8 x 10" 10 


-Mi/Mi 


0.32 - 2.0 x 10~ 3 


0.51 - 3.2 x 10" 2 


BH mass (M ) 


> 6 


> 8.6 



evolutio nary model, which is valid fo r M c < O.45M , has been ap- 
plied bv lKind dl993: hereafter K93h to thecase of the BH LMXB 
V404 Cyg with convincing results. Using the orbital period of GX 
339-4 (P = 1.7557 d) into the K93's equations we obtain 

(m c /0.25) 7 ' 65 m2 5 = 0.107 (5) 

with m c = M c /Mq and m 2 = M 2 /Mq. According to K93 M c 
is constrained between 

M 2 >M C > 0.17M 2 

where the right hand side condition comes from the Schonberg- 
Chandrasekhar limit. 



3.1 Minimum and maximum mass for the companion 

If we use the condition m c = m2 in eq(5]we obtain M c = M 2 = 
O.166Af . This mass is the minimum mass permitted for the com- 
panion and represents the case in which all the Hydrogen shell is 
burned and the companion becomes a Helium white dwarf. As we 
show in fig[T]as the dashed line, this limit case constrains the mass 
of the BH in GX339 to Mx > 6M , although this lower limit 
could be as high as Mx > 7.2Af if we consider the H03 upper 
limit to f(M) (i. e. 6.3M ) in the K-correction computation. 
On the other hand, we obtain M 2 < 1.1M (Af c < O.187M ) 
by applying the condition M c /M 2 ~ 0.17. This is the maximum 
mass permitted for the companion through the stripped-giant model 
and sets an upper limit to the mass ratio of q < 0.125 (fig. [5}. If the 
companion were as massive as 1.1M the BH would be more mas- 
sive than 8.6M . We list the solutions of the stripped-giant model 
for GX339-4 in TableQ] 



4 OBSERVATIONAL CONSTRAINTS TO GX 339-4 

Following K93, the radius, luminosity and mass transfer rate of a 
stripped-giant Roche lobe filling companion can be calculated as 
function of the mass of the Helium core according to: 

R 2 = 12.55i? (m c /O.25) 5 1 (6) 

L 2 =33L (m c /O.25) 811 (7) 

- Mi = 5.4 x 10 _9 (m c /0.25) 7 "ma (8) 

For the minimum mass solution we obtain R 2 = 1.56f? and 
L 2 — 1.19L which gives a companion effective temperature 
T e ff = 4837 K. On the other hand, if we consider the maxi- 
mum mass solution we obtain R 2 = 2.93i? , L 2 = 3.25L and 



K-correction for GX 339-4 3 



Table 2. Companion star in XTE J1550-564 





Orosz et al. 2002 


Stripped giant model 


Mass (Mq) 


0.4-1.3 


0.16-1.09 


Radius (Rq) 


1.88-2.81 


1.42-2.67 


Luminosity (Lq) 


1.4-4 


1 -2.8 


Teff(K) 


4100-5100 


4574-4880 


-M 2 (Mq yr- 1 ) 


0.1 - 7 x 10 -9 


0.04 - 0.7 x 10" 9 



T e f f = 4533 K. Although spectroscopic features of the compan- 
ion star have not been detected yet, SFC01 established a lower limit 
to the r-band Gunn magnitude of the companion of r > 20.4. We 
have computed the expected r magnitudes for our derived com- 
panions as function of the distance (d) to GX 339-4 which is es- 
timated to be in the range 6 < d < 15 kpc (H04). For the 
calculation we have used the mean value for th e reddening of 
E(B - V) ~ 1.2 reported in H04. According to lSchmidt-Kaled 
dl982[) , Ro = R + 2.32E(B - V) where R is the dereddened 
magnitude which is computed foll owing the r to R transformation 
reported in Windh orst etafl jl99lh . For d ~ 6 kpc we estimated 
21.8 > r > 20.6 for the minimum and maximum mass solutions 
respectively whereas we derive 23.8 > r > 22.6 if d ~ 15 kpc. 
Therefore, it is clear that a stripped-giant companion is compatible 
with the non-detection of the secondary by SFC01. 
On the other hand, eq[8] gives a mass transfer rate of 4.9 x 
10 _11 A/ e yr -1 and 7.8 x 1O" 1O M yr _1 for the minimum and 
maximum mass solutions respectively. KKB96 showed that the 
critical mass transfer rate (M c t) for a system to be a persistent 
source (i. e. not transient) can be estimated by using the expres- 
sion 

M ct = 5 x l O - 11 (m 1 ) 2/3 (-^) 4/3 M yr- 1 (9) 
ihr 

where mi is the mass of the compact object in M and the or- 
bital period in units of 3 hours. If we substitute in this equation the 
orbital period of GX 339-4 and assuming Mx ~ lOAf we obtain 
Met ~ 8 x 1O _9 M yr -1 , which shows that our values obtained 
through the stripped-giant model are consistent with the transient 
nature of the s ource . 

Homa net alj J2005t) detected a maximum X-ray flux during the 
2002/2003 outburst of Fx = 3.25 x 10 -8 erg s _1 cm" 2 which 
results in a X-ray luminosity in the range Lx (peak) ~ 1.4 — 8.7 x 
10 38 erg s _1 depending on the assumed distance. Applying 

L x (peak) = T)hhc (10) 

withr; ~ 0.1 jFrank, King. & Rainell992h we obtain M ~ 0.25- 
1.54 x 1O _7 M yr" 1 . Therefore, we find a predicted outburst duty 
cycle for the outburst in the range -M 2 /Mi ~ 0.32 - 2.0 x 10" 3 
for the minimum mass solution and —M 2 /M\_ ~ 0.51 — 3.2x 10 -2 
for the maximum companion mass case. Classical transient sources 
(e. g. V404 Cyg) only show an outburst every few decades and 
hence its duty cycles are low (~ 10~ 4 ). This is clearly not the be- 
haviour showed by GX 339-4 which has undergone four outbursts 
in the last ~ 10 years. This suggests that M 2 is close to M c t for 
GX 339-4, and therefore a companion mass in the upper part of the 
proposed mass range seems to better explain the high X-ray activity 
of this system. 
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5 XTE J1550-564, A TWIN SYSTEM OF GX 339-4? 

The LMXB XTE J1550-564 was discovered on 1998 by the All 
Sky Monitor (ASM) onboard the Rossi X-ray Tim ing Explorer 
JSmith et alJI 19981) . lOrosz et al.1 J2002; hereafter O02h showed that 
this system harbours a BH with f(M) = 6.86 ± 0.71M Q and 
Mx ~ IOMq. These authors also classify the companion star 
in the range G8IV-K4ni with 4100 < T eff < 5100 K and 
M 2 ~ O.4M0 although a tentative measurement of the rotational 
broadening suggests M 2 > 0.9A/q. On the other hand, the orbital 
period is P = 1.552 days (O02) which is very similar to that of 
GX 339-4. The stripped-giant model (see above equations) for this 
P yields 0.16 < M 2 < 1.09M Q with 4880 > T eff > 4574 
K, which is highly consistent with the observations. In table [2] we 
compare the observed properties of the companion in XTE J1550- 
564 with those obtained through the stripped-giant model and the 
parameters are in excellent agreement. 

Since its discovery in 1998, XTE J1550-564 has undergone two 
more X-ray outbursts suggesting a mass transfer rate close to M c t . 
O02 estimate M 2 = 0.1 - 7 x 1O~ 9 M yr" 1 which yields 
M 2 /M ct ~ 1.5 x 10~ 2 - 1. During the 1998 outburst. iHoman et all 
d200lh measured a maximum X-ray flux Fx = 1-47 x 10 erg 
s" 1 cm -2 which results in Lx(peak) ~ 4.94 x 10 38 erg s _1 
(d ~ 5.3 kpc; O02) and Mi ~ 8.72 x 10~ 8 A/ Q yr -1 (eqUOjl. By 
combining the observational constraints to M 2 and Af 1 we obtain a 
predicted duty cycle in the range -M2/M1 ~ 0.12-8.03 x 10" 2 . 
As expected, these values are larger that those typical of classi- 
cal transients (e. g. V404 Cyg) which have longer outburst re- 
currence times. We get to reproduce these values by considering 
stripped-giant companions with M 2 > 0.28, for which we pre- 
dict —A/2 in agreement with the observations. We want to note 
that XTE J1550-564 is also interesting because is one of only four 
SXTs for which the Very High X-ray state has been observed (e. g. 
Sobezak et al.1 19991; GX 339-4 is one of the other three). Moreover 
Hannikaine n et alj d200ll) observed relativistic plasma ejections at 
radio wavelenghts indicating that XTE J 1550-564 is a microquasar 
as well. 



6 DISCUSSION 

We have applied the stripped-giant model to the BH LMXBs GX 
339-4 and XTE J 1550-564. These systems share the following 
properties: 

(i) Orbital periods in the range 1.5-1.7 d. 

(ii) Transient behaviour with frequent X-ray outburst 

(iii) The Very High state has been observed (a total of only 4 
system have shown this X-ray state) 

(iv) Both are microquasars. 

For the case of XTE J1550-564, where the companion star has 
been detected and its stellar parameters constrained, we find that 
the proposed evolutionary model successfully reproduces all the 
observables. Therefore, it seems probable that XTE J1550-564 
harbours a stripped-giant companion. The case of GX339 is 
less straightforward since the companion star has not yet been 
observed. Only through the detection of the NIII/CIII Bowen lines 
arising from the irradiated donor a lower limit (K em ~ 317 km 
s _1 ) to the velocity of the companion was established by H03. We 
have applied the K-correction to this K em velocity and derived 
Mx > 6Mq including the error bars in /(M) reported by H03. 
This value represents a solid lower limit to the mass of the BH 
in this LMXB. This result comes from the assumption that the 



companion is a stripped-giant with the minimum possible mass. 
Although there is not a definitive evidence for a stripped-giant 
donor in this system, this model is favoured by the ~ 1.7 d orbital 
period which clearly points to an evolved companion. Moreover 
this evolutionary model is consistent with the transient behaviour 
of the source and also with the large number of X-ray outbursts dis- 
played. The radius and luminosity predicted by the stripped-giant 
model also explain the non-detection of the companion by SFC01. 
However, the above lower limit to Mx is quite conservative and 
unrealistically small since the companion is not an Helium white 
dwarf because it would not fill a 1.7 d Roche lobe. Furthermore, 
the M 2 predicted by the minimum mass solution yields an outburst 
duty cycle which is much too long to explain the frequent X-ray ac- 
tivity in GX 339-4, with 4 outbursts in the last 10 years. Assuming 
the ratio M 2 /M ct for GX 339-4 is similar than for XTE J1550-564 
{Mil M ct > 1-5 x 10~ 2 ) we estimate M 2 > 1.2 x 10" 10 M o 
yr -1 , which results in M 2 > 0.3A/q. If we combine this limit 
with the K-correction we obtain Mx > 6.6 — 7.7M© considering 
the error bar in f(M). Although this lower limit is less secure than 
the one obtained through the minimum mass solution it is probably 
more realistic and better explains the frequent X-ray activity 
displayed by this source. Note that the stripped-giant model is 
only valid for M c < 0.45M© and hence if we consider heavier 
evolved companions (e. g. giant stars) with M 2 > M c >O.45M0 
we would obtain higher Mx values. However, as we explain in 
section [3] this possibility is at odds with the mean density of the 
donor derived from the orbital period. 



7 CONCLUSION 

We have applied the K-correction to the BH LMXB GX 339-4. 
By considering the limit case where the emission line is formed 
at the limb of the irradiated region of the companion we derive 
a solid lower limit to Mx > 6M0 including the error bars 
in f(M). Here we have only assumed M 2 > O.I66M0, the 
lower limit allowed by the stripped giant model. We find that the 
stripped-giant evolutionary model explains the non-detection of 
the companion by SFC01 and the X-ray behaviour of the source. 
In particular, we propose M 2 > O.3M0, for which we predict M 2 
large enough to explain the frequent X-ray outbursts displayed by 
this source. This limit results in Mx > 7Mq . From the maximum 
mass solution we find q < 0.125. On the other hand, we have 
also shown that the stripped-giant model successfully explains the 
observable properties (M 2 , R 2 , T e ff and M 2 ) of the akin LMXB 
XTE J 1550-564. 
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